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ABSTRACT 

We present FUSE observations of the hot white dwarfs in the post-common envelope binaries 
Feige 24, EUVE J0720-317, BPM 6502, and EUVE J2013+400. The spectra show numerous pho- 
tospheric absorption hues which trace the white dwarf orbital motion. We report the detection of 
C III, O VI, P V, and Si IV in the spectra of Feige 24, EUVE J0720-317 and EUVE J2013-)-400, 
and the detection of C III, N II, Si III, Si IV, and Fe III in the spectra of BPM 6502. Abundance 
measurements support the possibility that white dwarfs in post-common envelope binaries accrete 
material from the secondary star wind. The FUSE observations of BPM 6502 and EUVE J2013+400 
cover a complete binary orbit. We used the FUSE spectra to measure the radial velocities traced 
by the white dwarf in the four binaries, where the zero-point velocity were fixed using the ISM 
velocities in the line of sight of the stellar systems. For BPM 6502 we determined a white dwarf 
velocity semi-amplitude of -ftTwD = 18.6 ± 0.5km s~^, and with the velocity semi-amplitude of the 
red dwarf companion (-ftTRD = 75.2 ± 3.1km s^^), we estimate the mass ratio to be g = 0.25 ± 0.01. 
Adopting a spectroscopic mass determination for the white dwarf, we infer a low secondary mass of 
Mrd = 0.14 ± 0.01 Mq. For EUVE J2013-I-400 we determine a white dwarf velocity semi-amplitude 
of KvfY) = 36.7 ± 0.7km s'^. The FUSE observations of EUVE J0720-317 cover approximately 30% 
of the binary period and combined with the HST GHRS measurements (Vennes et al. 1999, ApJ 523, 
386), we update the binary properties. FUSE observations of Feige 24 cover approximately 60% of 
the orbit and we combine this data set with HST STIS (Vennes et al. 2000, ApJ, 544, 423) data to 
update the binary properties. 
Subject headings: stars: abundances — binaries: spectroscopic — white dwarfs 



1. INTRODUCTION 

Post common-envelope binaries consist of an evolved 
primary (white dwarf or sdB) and a late-type main- 
sequence secondary in close orbit. These binaries are 
thought to have evolved from wide binary systems, where 
the more massive star evolved off the main sequence fill- 
ing its Roche lobe and beginning mass transfer onto its 
less massive companion. If the transfer is dynamically 
unstable, a common-envelope (CE) is formed and friction 
between the stellar components and the CE decreases the 
orbital separation and induces the ejection of the CE. 
Depending on the separation of the components, some 
of these systems further evolve to become cataclysmic 
variables. From the sample of well studied post-CE bi- 
naries (Schreiber & Gansicke 2003), approximately half 
will evolve into cataclysmic variables within a Hubble 
time. 

The atmosphere of a white dwarf a in close binary sys- 
tem usually displays enhanced traces of heavy elements. 
Vennes et al. (1999) found that the large abundance of 
carbon and helium as well as time variable helium abun- 
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dance in some white dwarfs provide evidence for on-going 
accretion from the red dwarf. 

We present high-resolution far-ultraviolet (FUV) spec- 
troscopic observations of Feige 24 (PC 0232+035: 
Vennes & Thorstensen 1994b; Vennes et al. 2000), 
WD 0718-316 (EUVE J0720-317, 2RE J0720-318: 
Vennes & Thorstensen 1994a; Barstow et al. 1995b), 
WD 1042-690 (BPM 6502: Kawka et al. 2000) and 
WD 2011-^398 (EUVE J2013-h400, 2RE J2013+400: 
Thorstensen et al. 1994; Barstow et al. 1995a) in §2. We 
present the analyses of the FUSE spectra, in §3, determine 
new orbital and stellar parameters of the four binary sys- 
tems in §4 and measure the heavy element abundance in 
the white dwarf atmospheres in §5. We summarize in §6. 

2. OBSERVATIONS 

We have obtained high-resolution FUV spectra of four 
close binary systems with FUSE (Table 1). The spectro- 
graph covers the FUV spectral range from 905 to 1187 A 
with a spectral resolution R = 20, 000±2000. The instru- 
ments are described in detail by Moos et al. (2000) and 
Sahnow et al. (2000a, b). The observations were made us- 
ing the LWRS (J0720-317, BPM 6502 and J2013-H400) 
and MDRS (Feige 24, J2013-H400) apertures and in time- 
tagged mode, except for Feige 24 which was observed in 
the HIST mode. The data for J0720-317 and BPM 6502 
were processed with the CALFUSE pipeline v3.1 and the 
data for Feige 24 and J2013-I-400 were processed with 
the CALFUSE pipeline v3.0. For the abundance anal- 
ysis, we co-added the individual exposures after align- 
ing them on the photospheric lines using the calculated 
radial velocities. Figure 1 shows the FUSE spectra of 
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Feige 24, J0720-317, BPM 6502 and J2013+400 indi- 
cating key heavy elements. Table 2 lists important pho- 
tospheric lines in Feige 24, J0720-317, and J2013+400, 
while Table 3 lists important lines in BPM 6502. 

3. WHITE DWARF ATMOSPHERIC PARAMETERS 

The FUSE spectrum of Feige 24 was analyzed by 
Vennes et al. (2005). Using LTE models they deter- 
mined an effective temperature of 64700 ± 3000 K and 
log.g = 7.58 ±0.25. Vennes et al. (2005) also investigated 
the effect of using LTE models on temperatures and sur- 
face gravities obtained from the analysis of Lyman and 
Balmer lines of hot white dwarfs. They found that for 
objects hotter 50000 K and assuming low metallicity, the 
LTE determinations are overestimated by approximately 
4000 K. Using the NLTE correction vector, the tempera- 
ture would be revised to approximately 60 000 K. Which 
within uncertainties is in agreement with previous tem- 
perature determinations, and in this work we will assume 
an effective temperature of Tea = 57 000 ±2000 K, which 
is based on estimates from previous spectroscopic studies 
(e.g., Vennes & Lanz 2001). 

The FUSE spectrum of BPM 6502 displays Lyman /3 
and 7 satellites. Hebrard et al. (2003) compared the 
FUSE spectrum to a LTE model spectrum at Tcs — 
21 380 and logg = 7.86 that included the quasi-molecular 
satellites of Lyman a, /3 and 7. The theoretical spectrum 
showed a reasonably good agreement with the observed 
FUSE spectrum of BPM 6502, however some discrepan- 
cies are observed. Hebrard et al. (2003) noted that their 
models do not include the variation of the dipole mo- 
ment during the collision, which may have a significant 
effect on the strengths of the satellite profiles. In this 
work, we will adopt the effective temperature and sur- 
face gravity determined by Kawka et al. (2007). Even 
though these parameters were determined using optical 
spectra, only Ha is significantly contaminated, and as a 
result was excluded in their analysis. 

We used the FUSE spectra to obtain an effective tem- 
perature, surface gravity and helium abundance of the 
white dwarfs in J0720-317 and J2013+400. We com- 
puted a grid of LTE plane-parallel models. The grid 
of models extends from T^s = 30 000 to 70 000 K (in 
steps of 4000 K), logg = 7.0 to 9.5 (in steps of 0.25 dex) 
and log(iVHo/A^H) = ~4.0 to 0.0 (in steps of 0.5 dex). 
We fitted 6 channels (SICIB, SIC2A, LIF2B, LIFIA, 
SICIA and SIC2B) simultaneously using a x^ minimiza- 
tion technique to obtain an effective temperature, sur- 
face gravity and helium abundance. Regions which show 
interstellar absorption features were excluded from the 
fit. We obtained T ^ 52 750 ± 150, logg = 7.73 ± 0.03 
and log (A^He/iVn) = -3.28 ± 0.08 for J0720-317. For 
J2013+400 we obtained T = 47 800 ± 200, logg = 
8.20 ± 0.03 and log {NujNn) = -2.90 ± 0.08. The 
atmospheric parameters derived from FUSE spectra for 
J0720— 317 are consistent with parameters derived from 
Balmer line spectroscopic fits. The effective tempera- 
ture and helium abundance for J20 13+400 are consis- 
tent with optical spectral analyses, however the surface 
gravity is significantly higher than the values determined 
from Balmer line spectral fits. This higher surface grav- 
ity, which corresponds to a mass of 0.80±0.02 Mq (using 
the mass-radius relations of Wood (1995)) is also signifi- 
cantly higher than the mass (0.64 ±0.03 M©) determined 



from the gravitational rcdshift of Vennes et al. (1999). 
We will adopt the optical values in this work until this 
discrepancy is resolved. We discuss this problem further 
in §4.4. 

Table 4 summarizes the atmospheric properties of the 
white dwarfs determined from FUSE spectra. The Table 
also compares these white dwarf parameters to the same 
parameters based on optical spectroscopic studies. 

4. BINARY PARAMETERS 

We used the FUSE spectra to obtain radial velocities 
of the four white dwarfs in the binary systems, Feige 24, 
J0720-317, BPM 6502 and J2013+400. These new ra- 
dial velocity measurements are given in Table 5. We 
phased the white dwarf radial velocities using published 
binary ephemerides and determined the velocity semi- 
amplitude K-ifi/D and mean velocity 7wd- These new 
parameters are then combined with the results of pre- 
vious optical studies of the red dwarf companions which 
listed the red dwarf velocity semi-amplitudes (-RTrd) and 
mean velocities (7rd)- From this set of measurements 
we calculated the binary mass ratios and the white dwarf 
gravitational redshifts (wg^wo) for all four systems. Fig- 
ure 2 shows the radial velocity measurements of the four 
binaries folded on their orbital period. 

Two small corrections were applied to the red dwarf 
data. First, a determination of the systemic velocity 7sys 
requires that we subtract an estimate of the gravitational 
redshift of the red dwarf (wj^rd) from its mean velocity. 
Secondly, a correction is applied to the red dwarf veloc- 
ity semi-amplitude to take into account the non-uniform 
distribution of the Ha emission over the red dwarf sur- 
face. Both corrections require initial estimates of the 
red dwarf mass and radius. These estimates could be 
obtained by adopting a mass and radius corresponding 
to published companion spectral type, but we favor the 
following approach. 

The strength of the Ha emission line has been observed 
to vary in all four systems (Vennes & Thorstensen 1994b; 
Kawka et al. 2002; Thorstensen et al. 1994) and is max- 
imum near superior conjunction. These cyclical varia- 
tions are explained by the changing viewing angle of the 
irradiated red dwarf hemisphere over the binary period 
(Thorstensen et al. 1978). In these systems, the Ha emis- 
sion originates from the irradiated hemisphere and traces 
an orbit that is lower than the true center of mass orbit. 
The irradiating flux originates from the hot photosphere 
of the white dwarf (e.g., Feige 24), or, alternatively, from 
a hypothetical hot accretion region on the surface of the 
white dwarf as proposed by Maxted et al. (1998) in the 
case of GD 448. 

We correct for the irradiation effect using the formal- 
ism described by Vennes et al. (1999). Wade & Home 
(1988) and Orosz et al. (1999) discuss a calculation as- 
suming an irradiated hemisphere that employs the same 
formalism. The correction term uses a preliminary esti- 
mate of the red dwarf radius obtained using a spectro- 
scopic determination of the white dwarf mass and the 
uncorrected mass ratio. The red dwarf radius is then 
calculated using the mass radius relations of Caillault & 
Patterson (1990). The procedure may be iterated at will. 
However, the initial correction to the mass ratio being of 
the order of a few percent (5 to 8%), further iterations 
are futile. 
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FUSE Observation Log 






Star 


Date 


Exposure time 


Aperture 


Dataset 


Observer 


Feige 24 


2003 Dec 07 


2149 


MORS 


P1040501 


Moos 




2004 Jan 02 


16117 


MORS 


P1040503 


Moos 




2004 Jan 06 


11999 


MORS 


P1040504 


Moos 


J0720-317 


2001 Nov 13 


17700 


LWRS 


B0510101 


Vennes 


BPM 6502 


2002 Jun 26 


32100 


LWRS 


Z9104501 


Andersson 




2006 Jan 21 


10696 


LWRS 


U1074601 


Blair 




2006 Jan 22 


7122 


LWRS 


U1074602 


Blair 




2007 Mar 02 


12758 


LWRS 


U1074604 


Blair 


J2013+400 


2000 Oct 11 


11200 


LWRS 


P2040401 


Moos 




2002 Oct 29 


6543 


LWRS 


M1053101 


Dupuis 




2003 Oct 21 


42607 


MORS 


D0580101 


Vidal-Madjar 




2003 Oct 23 


12881 


LWRS 


M1053102 


Dupuis 
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FlG. 1.— FUSE spectra of Feige 24, J0720-317, BPM 6502 and J2013+400 showing the key heavy elements and ISM lines. 
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Radial velocities of the white dwarf in Feige 24, J0720-317, BPM 6502 and J2013+400 folded on the orbital period as described 



Kawka et al. 



TABLE 2 

Photospheric lines in Feige 24, 

J0720-317, AND J2013+400. 



Ion A (A) 



Ion A (A) 



N IV 921.994 
N IV 922.519 
N IV 923.676 
N IV 924.284 
S VI 933.378 
S VI 944.523 
P IV 950.657'' 
N IV 955.334 
C III 977.020 
P IV 1030.515=^ 
O VI 1031.912 
P IV 1033.112=^ 
P IV 1035.516=^ 
O VI 1037.613 
S IV 1062.664 
Si IV 1066.63'= 
S IV 1072.973 
S IV 1073.518 



He II 1084.94<: 
Si III 1108.358 
Si III 1109.97^^ 
Si III 1113.23'' 
P V 1117.978 
P IV 1118.552" 
S V 1122.042 
Si IV 1122.485 
P V 1128.008 
Si IV 1128.34^ 
S V 1128.667 
C III 1174.993 
C III 1175.263 
C III 1175.590 
C III 1175.711 
C III 1175.987 
C III 1176.370 



''In Feige 24 only. P IV 950.657 may 
be blended with ISM O I 950.885 

•"Blend of 1066.614, 1066.636 and 
1066.650 A. 

•^In J0720-317 and J2013+400 only. 

<*Blend of 1109.940 and 1109.970 A. 

''Blend of 1113.174, 1113.204 and 
1113.230 A. 

^Blcnd of 1128.325 and 1128.340 A. 



TABLE 3 

Photospheric lines in 

BPM 6502. 



Ion A (A) 



Ion A (A) 



N II 1083.990 
N II 1084.576=^ 
N II 1085.545*= 
N II 1085.701 
Si III 1108.358 
Si III 1109.97^= 
Si III 1113.23^* 
Si IV 1122.485 
Fe III 1122.526 
Fe III 1124.881 
Fe III 1126.729 



Fc III 1128.050 
Si IV 1128.34° 
Fe III 1128.724 
Fe III 1129.191 
C III 1174.933 
C III 1175.263 
C III 1175.590 
C III 1175.711 
C III 1175.987 
C III 1176.370 



''Blend of 1084.562 and 1084.580 

A. 

•"Blend of 1085.529 and 1085.546 



'^Blcnd of 1109.940 and 1109.970 



A. 

A. 

^Blend of 1113.174, 1113.204 and 
1113.230 A. 
■"Blend of 1128.325 and 1128.340 

A. 



Final estimates of the red dwarf mass and radius 
are calculated from the corrected mass ratio and the 
mass of the white dwarf and from the mass-radius rela- 
tions of Caillault & Patterson (1990). Next, the grav- 
itational redshift of the red dwarf is calculated using 
Vg = O.636O8M(M0)/i?(i?0) km s^^ and is typically of 
the order of 0.5 km s^^. Finally, the systemic velocity 
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Fig. 3. — 2MASS J — H versus H — Ks color diagram showing 
the colors of the secondary stars (corrected for the white dwarf 
contribution) in Feige 24, J0720-317, J2013-I-400 and BPM 6502 
{open squares), compared to main-sequence colors (filled circles). 



is calculated using 7sys = 7rd — i'g,RD- Recent measure- 
ments of late-type main sequence stars using interferom- 
etry and eclipsing binaries show a relatively large scat- 
ter in the mass-radius relations for these stars (Berger 
et al. 2006; Lopez- Morales 2007). Rebassa-Mansergas 
et al. (2007) show that a similar scatter is observed in 
the late-type main sequence companions to white dwarf 
stars. However, the effect of this scatter on the gravita- 
tional redshift of the secondary is less than 0.1 km s~^ 
and does not affect our results. 

We estimate the spectral type of the secondary stars 
using the mass- type relations of Kirkpatrick & McCarthy 
(1994), but we corroborate our spectral type determina- 
tion using 2MASS JHKg photometry. First, we calcu- 
late the absolute JHK magnitudes of the white dwarfs. 
These white dwarf magnitudes are then converted from 
CIT to 2MASS using Cutri et al. (2006)^ Next, we cal- 
culate the absolute magnitude of the system using the 
geometric parallax (Feige 24) or an ultraviolet based pho- 
tometric parallax of the white dwarf uncontaminated by 
the companion (J0720-317, BPM 6502 and J2013-h400). 
Next, we subtract the white dwarf contribution. Figure 3 
shows the 2MASS JHKg of the secondary companions 
in the binary systems, Feige 24, J0720-317, BPM 6502 
and J2013-I-400, compared to the colors for cool main- 
sequence stars (Bessell & Brett 1988; Bessell 1991) con- 
verted to 2MASS using Cutri et al. (2006). 

Finally, we determine the gravitational redshift of the 
white dwarf (wg) by subtracting the systemic velocity 
(7sys) from the white dwarf mean velocity (7wd)- Then, 
we convert the gravitational redshift of the white dwarf 
into a mass estimate using mass-radius relations of Wood 
(1995). 

With this general frame of work in mind, we now dis- 
cuss each system separately. The binary properties of 
these systems are given in Table 6, while the properties 
of the component stars are given in Table 4. 

® Available at 
http://www.ipac.caltcch.edu/2mass/releases/allsky/doc/cxplsup.html 
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TABLE 4 
Properties of the white dwarf and its cool companion. 



Parameter'' 



Feige 24 



J0720-317 



BPM 6502 



J2013+400 



White Dwarf 



Tcff.opt (K) 


57 000 ±2000'' 


52 400 ± 1800<= 


19 960 ± 400^* 


48 000 ± 900<= 


Teff.FUSE (K) 


f 


52 750 ± 150 


... g 


47 800 ± 200 


logSopt (c.g.s.) 


7.66 ± 0.08 


7.68 ±0.01== 


7.86 ± 0.09^ 


7.69 ± 0.09° 


logflFUSE (c.g.s.) 


f 


7.73 ± 0.02 


... g 


8.20 ±0.03 


log{ArHc/AfH) 




-3.28 ±0.08 




-2.90 ±0.08 


Mopt (Mq) 


0.58 ±0.05^ 


0.56 ± 0.04== 


0.55 ±0.05 


0.56 ± 0.03° 


Vg (km s"-') 


20.1 ±1.9 


21.4 ±1.9 


17.9 ±0.5 


34.0 ±1.3 


MgR (Mq) 


0.57 ±0.03 


0.58 ±0.03 


0.46 ±0.01 


0.71 ±0.02 


Red Dwarf 


M {Mq) 


0.39 ± 0.02 


0.43 ± 0.03 


0.14 ±0.01 


0.23 ±0.01 


K{Rq) 


0.43 ± 0.02 


0.47 ±0.03 


0.19 ±0.02 


0.29 ±0.01 


Vg (km s ^) 


0.6 ±0.1 


0.6 ±0.1 


0.5 ±0.1 


0.5 ±0.1 


Spec. Type. 


M2 


M2 


M5 


M3 


^AU parameters are 


determined in this work unless indicated otherwise. 




''From Vennes & Lanz (2001). 








^From Vennes et al. 


(1997). 








From Kawka et al. 


(2007). 








'^From Vennes ct al. 


(1999). 









See text. 
^See Hcbrard et al. (2003) for a comparison between optical parameters and FUSE spec- 
trophotometry (see text for details). 

Determined using the measured parallax. 



4.1. Feige 24 

Thorstensen et al. (1978) were the first to obtain or- 
bital parameters for Feige 24 by measuring the red dwarf 
radial velocities using Ha and He I emission lines. Vennes 
& Thorstensen (1994b) determined an orbital period of 
4.23160 ± 0.00002 days and using the absorption lines 
they measured a red dwarf velocity semi-amplitude of 
75.5 ± 2.1 km s~^ with a red dwarf mean velocity of 
7rd = 62.0 ± 1.4 km s~^ . Therefore, the measurements 
represent the true orbit traced by the secondary center 
of mass. 

Vennes et al. (2000) obtained HST STIS spectra at or- 
bital quadratures and obtained an estimate of the white 
dwarf semi-amplitude of 49.1 ± 0.3 km s^^ and a white 
dwarf mean velocity of 79.6 ± 2.3 km s~^. We measured 
the radial velocities of the white dwarf using the Si IV 
A1066.63 A line in the FUSE spectra. We used the ISM 
line of O I A1039.230 A to fix the zero-point of the wave- 
length calibration. We fixed the velocity of the ISM to 
+7.4 km s^^, which is the mean velocity of the two ISM 
components in the line of sight of Feige 24 (Vennes et 
al. 2000). We phased the FUSE data to the orbital 
ephemeris of Vennes & Thorstensen (1994b) to obtain 
a white dwarf semi-amplitude of 51.2 ± 0.6 km s^^ with 
a white dwarf mean velocity of 81.6 ± 0.4 km s"^. The 
FUSE and HST white dwarf mean velocities are in agree- 
ment within uncertainties, however the white dwarf semi- 
amplitude determined using FUSE is slightly larger than 
determined by HST. Note that the FUSE spectra lacks 
coverage at one of the quadratures ($ = 0.25). Combin- 
ing the FUSE a.nd HST data,, we obtain i^wD = 51.0±0.5 
km s~^ and 7wd — 81.5 ± 0.4 km s~^. Figure 2 shows 
the radial velocity measurements folded on the orbital 
period. 

A minimum white dwarf mass of 0.53 Mq is obtained 
using Kepler's third law. Vennes et al. (2000) estimated 



a white dwarf mass of 0.55 ± 0.02 Mq using the parallax 
measurement of tt = 14.7±0.6 mas (Benedict et al. 2000), 
V — 12.56 ± 0.05 for the white dwarf (Holberg et al. 
1986) and adopting Teff = 56 000 ± 1000 K. Using an 
effective temperature of 57 000 ±2000 K we obtain a mass 
of 0.58 ± 0.05 Mq. Note that the increased error in the 
temperature results in an increased error in the mass 
determination. 

We determined the white dwarf gravitational redshift 
to be Vg = 20.1 ± 1.9 km s~^. The gravitational redshift 
of the red dwarf was assumed to be 0.6 ±0.1 km s~^ (see 
below). The gravitational redshift of the white dwarf 
corresponds to a mass of 0.57 ± 0.03 M©. 

Since the two mass estimates agree, we adopt a mass 
of 0.57 ± 0.03 Mq for the white dwarf, which is the 
weighted mean of the two mass determinations discussed 
above. The corresponding surface gravity would be 
log 5 — 7.66 ± 0.08. The mass ratio of the system is 
q = i^wD/-f^RD — AfRD/AfwD = 0.68 which results 
in a mass of 0.39 ± 0.02 Mq for the secondary. And 
with a radius of 0.43 ± 0.02 Rq, we obtain a gravita- 
tional redshift of 0.6 ± 0.1 km s~^. The mass function of 
/(Mwd) — 0.189 ± 0.016 implies the inclination of the 
system is i = 77.4 ±3.5°. 

The derived mass (0.39 Mq) translates to a spectral 
type of dMl.5-2 for the secondary star. Using the 2MASS 
photometry of the binary system (J — 11.265 ± 0.024, 
H = 10.733 ±0.022, if, = 10.557 ±0.019), we calculated 
the absolute magnitudes of the secondary, Mj = 7.28, 
Mh = 6.65, and Mk — 6.45. In our calculations we 
used the parallax distance of 68.4 ± 2.0 pc (Benedict et 
al. 2000). Figure 3 shows that the spectral type of the 
secondary based on 2MASS colors {J — H = 0.63, H — 
Ks = 0.20) ranges from ~ K6-M2. 

The binary parameters and mass of the secondary 
(0.39 Mq) suggests that magnetic braking will be the 
main angular momentum loss for the system, and us- 
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TABLE 5 
FUSE Radial velocity measurements. 



HJD 

(2450000+) 



Phase 



(km s ^) 



HJD Phase v HJD 

(2450000+) (km s^^) (2450000+) 



Phase 



(km s ^) 









Fcige 24 










2980.556641 


0.306 


39.09 


3007.908447 


0.770 


134.65 


3010.681641 


0.425 


61.57 


2980.738525 


0.349 


40.54 


3007.925049 


0.774 


133.33 


3010.697021 


0.429 


56.18 


2980.743408 


0.350 


40.18 


3007.987061 


0.788 


136.60 


3010.749023 


0.441 


58.40 


2980.760498 


0.354 


43.99 


3008.000000 


0.792 


133.43 


3010.768555 


0.446 


67.56 


2980.765381 


0.355 


44.90 


3008.049072 


0.803 


130.67 


3010.820312 


0.458 


67.60 


3007.428223 


0.656 


122.77 


3008.063721 


0.807 


130.21 


3010.837158 


0.462 


68.66 


3007.443359 


0.660 


126.06 


3008.118164 


0.819 


128.69 


3010.905273 


0.478 


74.85 


3007.492432 


0.672 


127.64 


3008.133545 


0.823 


129.36 


3010.964844 


0.492 


76.92 


3007.509277 


0.676 


124.60 


3008.188477 


0.836 


123.89 


3010.980957 


0.496 


76.64 


3007.562744 


0.688 


125.24 


3008.271729 


0.856 


126.05 


3011.031006 


0.508 


78.30 


3007.578125 


0.692 


129.65 


3008.325439 


0.868 


119.41 


3011.044678 


0.511 


83.00 


3007.631592 


0.704 


127.55 


3008.341064 


0.872 


119.50 


3011.097412 


0.523 


90.51 


3007.706787 


0.722 


132.29 


3008.394775 


0.885 


116.56 


3011.118652 


0.528 


90.47 


3007.769775 


0.737 


133.51 


3008.410400 


0.889 


117.65 


3011.167725 


0.540 


92.46 


3007.786621 


0.741 


130.78 


3008.462891 


0.901 


109.57 


3011.183594 


0.544 


93.12 


3007.839355 


0.754 


133.54 


3010.547607 


0.393 


46.78 


3011.238281 


0.557 


97.22 


3007.855713 


0.757 


131.35 


3010.613037 


0.409 


54.23 


3011.251465 


0.560 


99.58 



J0720-317 



2226.887805 
2226.954556 



0.682 
0.735 



125.4 
129.5 



2227, 
2227, 



024896 
095490 



0.791 
0.847 



127.5 
112.2 



2227.164429 0.902 



93.7 









BPM 6502 










2452.329834 


0.742 


43.94 


2452.636963 


0.654 


40.07 


3758.507324 


0.115 


12.14 


2452.383789 


0.902 


37.55 


2452.708984 


0.868 


39.26 


4161.851074 


0.745 


46.36 


2452.438965 


0.066 


17.22 


3757.399658 


0.826 


42.38 


4161.914062 


0.933 


35.11 


2452.491455 


0.222 


6.87 


3757.468994 


0.032 


23.22 


4161.984863 


0.143 


13.60 


2452.562256 


0.432 


17.30 


3758.439453 


0.914 


35.99 


4162.055664 


0.353 


12.88 









J2013+400 










1858.616821 


0.261 


-31.31 


2934.230957 


0.833 


31.12 


2935.272705 


0.309 


-31.28 


1858.688843 


0.363 


-25.97 


2934.303467 


0.935 


16.50 


2935.344727 


0.412 


-18.91 


1858.759399 


0.463 


2.42 


2934.374268 


0.036 


-7.95 


2935.411865 


0.506 


12.85 


1858.827881 


0.560 


11.69 


2934.441895 


0.132 


-25.66 


2935.453857 


0.566 


21.00 


2576.795898 


0.206 


-28.09 


2934.483887 


0.191 


-32.40 


2935.482422 


0.606 


30.27 


2576.863037 


0.301 


-25.85 


2934.512695 


0.232 


-27.90 


2935.526123 


0.668 


33.93 


2576.932129 


0.399 


-9.26 


2934.557617 


0.296 


-32.68 


2935.553223 


0.707 


36.46 


2577.002197 


0.498 


12.66 


2934.631348 


0.400 


-15.26 


2935.599121 


0.772 


40.11 


2577.075684 


0.603 


34.30 


2934.712891 


0.516 


22.12 


2935.622559 


0.805 


35.33 


2936.028809 


0.381 


-21.45 


2934.767090 


0.593 


24.37 


2935.672607 


0.876 


28.87 


2936.099121 


0.481 


3.57 


2934.782471 


0.614 


29.15 


2935.689697 


0.900 


23.53 


2936.169678 


0.581 


29.15 


2934.834961 


0.689 


37.86 


2935.736572 


0.967 


14.54 


2936.241455 


0.682 


35.89 


2934.851807 


0.712 


40.96 


2935.805664 


0.065 


-16.10 


2936.312500 


0.783 


35.89 


2934.904785 


0.788 


40.11 


2935.822266 


0.088 


-21.72 


2936.383301 


0.884 


22.68 


2934.920898 


0.811 


36.46 


2935.875000 


0.163 


-25.94 


2934.023438 


0.539 


21.84 


2934.987793 


0.906 


23.53 


2935.891602 


0.187 


-31.84 


2934.082764 


0.623 


37.02 


2935.065674 


0.016 


-12.73 


2935.946045 


0.264 


-32.96 


2934.096191 


0.642 


34.77 


2935.125244 


0.100 


-23.97 


2935.960938 


0.285 


-32.40 


2934.156982 


0.728 


37.30 


2935.137695 


0.118 


-25.38 








2934.167969 


0.744 


40.11 


2935.199219 


0.205 


-32.96 









ing the equations from Schreiber & Gansicke (2003) and 
Rittcr (1986), wc find that the secondary star will fill 
its Roche lobe and begin mass transfer with a period of 
0.164 days in ~ 2.2 x 10^^ years. Therefore, Feige 24 is 
not representative of the progenitors of the current pop- 
ulation of cataclysmic variables. 

4.2. J0720~317 

Identified as a post-CE binary by Vennes & 
Thorstensen (1994a), Kawka et al. (2002) measured an 
orbital period for the binary of 1.262396 ± 0.000008 
days and a velocity semi-amplitude of the red dwarf of 
98.2 ± 1.2 km s~^ with a red dwarf mean velocity of 
7RD = 31.1 ±0.7 km s^i. Vennes et al. (1999) used HST 
spectra to trace the orbit of the white dwarf, and mea- 
sured the white dwarf semi-amplitude to be 79.5 ± 1.4 



km 



We have measured the radial velocities of the 



white dwarf using the Si IV A1066.63 fine in the FUSE 
spectra. We used the ISM line of O I A1039.230 A to 
fix the zero-point of the wavelength calibration. The 
velocity of the local interstellar cloud (LIC) in the di- 
rection of J0720— 317 is 13.5 km s~^ (Lallement et al. 
1995), however the measured velocity of O I A1039.230 
A is 6.7 km s~^. Therefore, all measured velocities were 
shifted by -1-6.8 km s""'^. Combining the FUSE velocity 
measurements with the HST velocities we updated the 
white dwarf semi-amplitude to ii'wD = 80.8 ± 1.2 km 
s^^ and 7wD = 51.9 ± 1.1 km s^^. Figure 2 shows the 
radial velocities of the white dwarf folded on the orbital 
period. 

We corrected the red dwarf semi-amplitude to obtain 
105.9 ± 3.4 km s~^, from which we estimate a mass ra- 
tio oi q = 0.76 ± 0.03. Using the white dwarf mass 
of 0.56 ± 0.04 Mq from Vennes et al. (1997b) and our 
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TABLE 6 
Binary Parameters 



Parameter 



Fcige 24 



J0720-317 



BPM 6502 



J2013+400 



my 

Period (days) 

/STrd (km s-l) 

^RD.corr (km S^^) 

TRD (km s^i) 
7sys (km s"^) 
A"wD (km s^i) 
7WD (km s^i) 
<? 

in 

/(mwo) {Mq) 



12.4 
4.23160 ± 0.00002=^ 

75.5±2.1'''d 

62.0 ±1.4=^ 

61.4 ±1.5 
51.0 ±0.5 

81.5 ±0.4 
0.68 ±0.02 

77.4 ±3.5 
0.189 ±0.016 



14.0 
1.262396 ±0.000008'' 
98.2 ±1.2*^ 
105.9 ±3.4 
31.1 ±0.7*= 
30.5 ±0.8 

80.8 ±1.2 

51.9 ±1.1 
0.76 ± 0.03 

72.2 ± 1.8 
0.155 ±0.015 



13.0 
0.3367849 ± 0.0000006 

71.1 ±0.2 

75.2 ±3.1 

8.7 ±0.1 

8.2 ±0.2 
18.6 ±0.5 
26.1 ±0.3 

0.25 ±0.01 

20.3 ±0.3 
0.015 ± 0.002 



14.8 

0.705517 ± 0.000006^= 

84.2 ± O.Q'^ 

89.1 ±2.6 

-29.5 ±0.7^^ 

-30.0 ±0.8 

36.7 ±0.7 

4.0 ±0.5 

0.41 ±0.01 

34.7 ±0.5 

0.052 ± 0.005 



*From Vennes & Thorstensen (1994b). 
''From Kawka et al. (2002). 
"^From Vonncs et al. (1999). 
Based on red dwarf absorption velocities. 



measured mass ratio, we estimate the mass for the red 
dwarf to be Mrd = 0.43 ± O.OSM© with a radius of 
0.47 ±0.03 i?0. Adopting the corrected mass function of 
/(Mwd) = 0.155 ± 0.015 Mq, we find the inchnation of 
the system to be i = 72.2 ± 1.8°. 

The derived mass translates to a spectral type of 
dMl-2 for the secondary star. Using the 2MASS pho- 
tometry of the binary system (J = 13.253 ± 0.025, 
H = 12.749 ± 0.026, Ks = 12.502 ± 0.027), we calcu- 
lated the absolute magnitudes of the secondary to be 
Mj = 6.95, Mh = 6.38 and Mr = 6.12. To esti- 
mate the distance toward J0720— 317, we obtained short- 
wavelength lUE spectra (SWP 54496,54497) and com- 
pared it to a synthetic spectrum with the white dwarf 
parameters (Tog- — 52400 K, log 5 = 7.68). The distance 
(d) to J0720— 317 would then be when the difference be- 
tween the observed spectrum and the model spectrum 
(F) placed at given distance, i.e., {R/d)'^F is minimized. 
The radius of the white dwarf (i?) is calculated using the 
mass-radius relations of Wood (1995). Figure 3 shows 
that the spectral type of the secondary based on 2MASS 
colors {J -H ^ 0.56, H - Ks = 0.26) is - dM3. 

We determined the white dwarf gravitational redshift 
to be Vg = 21.4 ± 1.9 km s^^. A red dwarf gravita- 
tional redshift of 0.6 ± 0.1 km s"'^ was used. The gravi- 
tational redshift of the white dwarf corresponds to a mass 
of 0.58 ± 0.03 Mq which is in agreement with the mass 
determined from a spectroscopic fit of the Balmer lines 
(Table 6). 

The binary parameters and the relatively large sec- 
ondary mass suggests that magnetic braking will be the 
main angular momentum loss for the system, and using 
the equations from Schreibcr & Gansicke (2003) and Rit- 
ter (1986), we conclude that the secondary star will fill 
its Roche lobe and begin mass transfer with a period of 
0.175 days (4.21 hrs) in - 3.2 x 10^ years. The con- 
tact period of 3.31 hrs reported by Vennes et al. (1999) 
is shorter than the one we determined in this work, be- 
cause of a larger secondary mass and hence higher mass 
ratio than the values used in Vennes et al. (1999). Due to 
the large mass ratio dynamically unstable mass transfer 
will be initiated when the system comes into contact (de 
Kool 1992). 
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Fig. 4. — Top: Periodogram of the radial velocities measured by 
Kawka et al. (2002). Middle: Periodogram of the radial velocities 
measured by Morales- Rueda et al. (2005). Bottom: Periodogram 
of the combined radial velocities. The lines shows the Icr (full line), 
2(7 (dashed) and 3a (dotted) confidence levels. 



4.3. BPM 6502 
4.3.1. Orbital Period 

Identified as post-CE binary by Kawka et al. (2000), 
Kawka et al. (2002) measured an orbital period of 
0.336784 ±0.000001 days, however Morales-Rueda et al. 
(2005) report a period of 0.337083 ± 0.000001 days that 
is significantly longer. The period determined by Kawka 
et al. (2002) corresponds to their third alias. We have 
reanalyzed these two data sets and we show that the cor- 
rect orbital period of BPM 6502 is the one reported by 
Kawka et al. (2002). 

The first step of our analysis was to determine the pe- 
riod of the two sets independently. The orbital period 
using measurements of Kawka et al. (2002) is 0.336784 ± 
0.000001 days, as reported. Aliases in the analysis are 
observed but all of them are significantly below the 3a 
confidence level. Figure 4 (top panel) shows the peri- 
odogram of these measurements and Figure 5 (top panel) 
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Fig. 5. — Radial velocities of the red dwarf in BPM6502 traced by 
Ha and folded on the orbital period as determined by Kawka et al. 
(2002) (top), Morales-Rueda et al. (2005) (middle) and the period 
determined from combining the radial velocity measurements from 
both studies (bottom). Velocities measured by Kawka et al. (2002) 
are indicated by open squares and velocities measured by Morales- 
Rueda et al. (2005) are indicated by full squares. 

shows Ha radial velocity measurements from both sets 
of data folded over the orbital period. Using the radial 
velocity measurements reported by Morales-Rueda et al. 
(2005), we reproduced the most probable orbital period 
of 0.337083 ± 0.000004 days, and the several aliases as 
reported in their analysis. We find that 4 possible peri- 
ods are within the 3a confidence level, including 2 that 
fall within la. The second peak that falls within la cor- 
responds to a period of 0.336786 days which corresponds 
to the period determined by Kawka et al. (2002) . There- 
fore, using the data of Morales-Rueda et al. (2005) alone, 
the remaining 3 aliases cannot be excluded as possible 
orbital periods of the system. Figure 4 (middle panel) 
shows the periodogram of these measurements, clearly 
showing that the other aliases cannot be excluded. The 
middle panel of Figure 5 shows Ha radial velocity mea- 
surements from both sets of data folded over the best or- 
bital period. Here the measurements taken during HJD 
24510143 - HJD 24510528 (Morales-Rueda et al. 2005) 
follow the calculated orbital velocities, however the mea- 
surements taken almost 3 years later HJD 2451607 - HJD 
2451634 (Kawka et al. 2000) and a further 2 years later 
HJD 2452301 - HJD 2452317 (Kawka et al. 2002) are out 
of phase. 

Using both sets of data we found that the best orbital 
period is: 

P = 0.3367849 ± 0.0000006 days 

and the epoch of inferior conjunction 

To = 2450143.4195 ±0.0003 

where the new ephemeris is formally in agreement with 
Kawka et al. (2002). Figure 4 (bottom panel) shows the 
periodogram of these measurements, where the aliases 
have almost disappeared. The bottom panel of Figure 5 
shows Ha radial velocity measurements from both sets 
of data folded over the best orbital period. The semi- 
amplitude of the combined measurements is -R'wd ~ 



71.1 ± 0.2 km s ^ with a red dwarf mean velocity, 
7RD = 8.7 ±0.1 kms-i. 

4.3.2. White Dwarf Radial Velocity Measurements 

We measured the radial velocities of the white dwarf 
using the narrow Si III lines Al 108.358, 1109.962, 
1113.225 A in FUSE spectra. We used the ISM hues 
of N I A1134.165, 1134.415, 1134.980 A to fix the zero- 
point of the wavelength calibration. All sets of veloc- 
ity measurements were corrected such that the ISM ve- 
locity is equal to uism = ^7 km s^^ (Lallcment et al. 
1995). The semi-amplitude of the white dwarf radial ve- 
locity measurements is Kwd = 18.6 ± 0.5 km s^^ and 
7WD = 26.1 ± 0.3 km s^^. Figure 2 shows the radial 
velocities of the white dwarf folded over the period. 

The Ha emission in BPM 6502 is assumed to have orig- 
inated from the irradiated hemisphere of the red dwarf 
(Kawka et al. 2002), and therefore we corrected the red 
dwarf semi-amplitude to 75.2 ± 3.1 km s^^. Using the 
corrected red dwarf semi-amplitude, we estimate a mass 
ratio oi q = 0.25 ± 0.01. The white dwarf mass of 
0.55 ± 0.05 Mq was estimated from Kawka et al. (2007) 
using the mass-radius relations of Wood (1995). From 
our measured mass ratio, we estimate the mass for the 
red dwarf to be Mrd — 0.14 ± O.OIM© and a radius of 
0.19 ±0.02 i?0. Adopting the corrected mass function of 
/(Mwd) = 0.015 ± 0.002 Mq, we find the inchnation of 
the system to be z = 20.3 ± 0.3°. 

The derived mass (0.14 Mq) translates to a spec- 
tral type of dM4.5(±0.5). Using 2MASS photometry 
of BPM 6502 (11.423 ± 0.026, H = 10.896 ± 0.027, 
Kg = 10.561 ± 0.021), we calculated the absolute mag- 
nitudes of the secondary, Mj = 8.88, Mh = 8.29 and 
Mk = 7.93. We determined the distance of BPM 6502 
using lUE spectra (SWP 27351) in the same way as for 
J0720— 317. Figure 3 shows that the spectral type of 
the secondary based on 2MASS colors [J — H = 0.60, 
H-K = 0.36) is ~ dM5. Tappert et al. (2007) obtained 
K-band spectroscopy, where the spectral energy distribu- 
tion suggests a spectral type ranging from M2.5 to M5, 
but the Na I and Ca I line strengths prefer a latter spec- 
tral type of M5. Taking into account all the available 
data we will adopt a spectral type of dM5. 

We determined the white dwarf gravitational rcdshift 



to be Vn 



17.9 ± 0.5 km s 



mass of 0.46 ± 0.01 M. 



which corresponds to a 

. A red dwarf gravitational red- 

^^ was used. This mass 



shift of Vg = 0.5 ± 0.1 km s 
is significantly different from the spectroscopic mass of 
0.55 ±0.05 Mq. Previous gravitational rcdshift measure- 
ments of the white dwarfs using Balmer lines resulted in 
much higher values, 42.2±9.0 km s~^ (Kawka et al. 2000) 
and 36.16±2.72kms-i (Morales-Rueda et al. 2005). It is 
likely that the Balmer line cores are contaminated by the 
red dwarf even after the subtraction of the emission line. 
In the case of using the FUSE spectra of the white dwarf 
to measure the radial velocity, the zero point depends 
on the adopted ISM velocity. We have assumed that the 
N I lines originate only from the LIC, however another 
cloud in the line of sight toward BPM 6502 could shift the 
velocity scale. Future high-dispersion ultraviolet obser- 
vations of the line of sight toward BPM 6502 may well be 
able to resolve velocity structures in ISM spectral lines 
and provide a more reliable anchor for ultraviolet line 
velocities. 
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The binary parameters and the relatively small sec- 
ondary mass suggests that major contributor to angular 
momentum loss is from the release of gravitational radi- 
ation. Using the equations from Ritter (1986), the sec- 
ondary star will fill its Roche lobe and begin mass trans- 
fer with a period of 0.083 days in ^ 3.0 x lO"'^" years. The 
time BPM 6502 comes into contact is much longer than 
for J0720— 317 because magnetic braking is not invoked 
due to the low-mass of the secondary. 

4.4. J 20 13+ 400 

Identified as post-CE binary by Thorstensen et al. 
(1994), Vennes et al. (1999) measured an orbital period 
for the binary of 0.705517 ± 0.000006 days and a semi- 
amplitude of the red dwarf of 84.2 ± 0.9 km s~^ with a 
mean velocity of 7rd = —29.5 ± 0.7 km s^^. Using HST 
spectra they were able to trace the orbit of the white 
dwarf, and measured the white dwarf semi-amplitude to 
be 29.9 ± 2.5 km s~^ . We have measured the radial ve- 
locities of the white dwarf using the Si IV A1066.63 line 
in the FUSE spectra. This line is near the ISM line of 
Ar I A1066.660, and blends with the Si IV A1066.63 line 
for a short duration near orbital conjunction. We only 
included the velocity measurements where the Si IV and 
Ar I are clearly separated. Ar I A1066.660 was used to 
fix the zero-point of the wavelength calibration. Its prox- 
imity minimizes the distortion in the wavelength scale 
due to thermal effects. All sets of velocity measurements 
were corrected such that the ISM velocity is equal to 
wisM = —8 km s~^ (Lallement et al. 1995). Combining 
the FUSE velocity measurements with the HST veloci- 
ties to obtain an updated white dwarf semi-amplitude of 
XwD ==36.7 ±0.7 km s-i and 7wd = 4.0 ± 0.5 kms^i. 
Figure 2 shows the radial velocities of the white dwarf 
folded over the orbital period. 

The Ha emission is assumed to have originated from 
the irradiated hemisphere of the red dwarf, and hence 
does not trace the motion of the red dwarf center of mass. 
Therefore, we corrected the red dwarf semi-amplitude to 
89.1 ± 2.6 km s""'^, from which we estimate a mass ra- 
tio oi q — 0.41 ± 0.01. Using the white dwarf mass 
of 0.56 ± 0.03 Mq from Vennes et al. (1999) and our 
measured mass ratio, we estimate the mass for the red 
dwarf to be Mrd = 0.23 ± O.OIM© and a radius of 
0.29 ±0.01 Rq. Adopting the corrected mass function of 
/(Mwd) = 0.052 ± 0.005 Mq, we find the inchnation of 
the system to be i = 34.7 ± 0.5°. 

The derived mass translates to a spectral type of 
dM3.5(±0.5) for the secondary star. Using 2MASS 
photometry of J2013+400 (J = 13.044 ± 0.044, H = 
12.520 ± 0.024, Ks = 12.260 ± 0.032), we calculated the 
absolute magnitudes of the secondary to be Mj = 7.73, 
Mr = 7.09 and Mk = 6.80. We determined the dis- 
tance of BPM 6502 using lUE spectra (SWF 27351) in 
the same way as for J0720— 317. Figure 3 shows that the 
spectral type of the secondary based on 2MASS colors 
[J -H = 0.64, H -Ks = 0.29) is ~ dM4. 

We determined the white dwarf gravitational redshift 
to be Vg — 34.0±1.3 km s^^, which corresponds to a mass 
of 0.71 ± 0.02 Mq. a red dwarf gravitational redshift of 
Vg ^ 0.5 ± 0.1 km s~^ was used. This is higher than the 
mass determined from a spectroscopic fit of the Balmer 
lines (see Table 6) and than the mass of 0.64 ± 0.03 Mq 
determined from the gravitational redshift by Vennes et 



al. (1999). Again, as in the case of BFM 6502, the pos- 
sibility that the adopted LIC velocity is blended with 
another ISM cloud causing a shift in the velocity scale 
could explain this difference. Moreover, the mass deter- 
mined from the FUSE spectral fit also indicates a higher 
mass (O.SOMq). The various mass estimates suggests 
that the white dwarf mass in J2013+400 is is between 
0.56 Mq and 0.80 Mq. 

The binary parameters and the relatively small sec- 
ondary mass suggests that major contributor to angu- 
lar momentum loss is from the release of gravitational 
radiation. Using the equations from Ritter (1986), the 
secondary star will fill its Roche lobe and begin mass 
transfer with a period of 0.118 days in ~ 1.3 x 10^^ years. 

5. WHITE DWARF ABUNDANCE ANALYSIS 
5.1. Model Atmospheres 

We computed a series of NLTE model atmospheres 
using TLUSTY version 200 and SYNSPEC version 48 
(Hubcny & Lanz 1995) and adopted stellar parameters 
for BFM 6502, J2013+400, J0720-317, and Feige 24: 
Teff = 20,000K (log .9 = 7.85), and T^s = 48,000 and 
52,400 K (logg = 7.7), and 57,000K (log^ = 7.5), re- 
spectively. Using preliminary abundance estimates, we 
then bracketed the abundance inputs with models with 
abundance variations of —0.7, and +0.7 dex. Table 7 
summarizes adopted model atoms, including some from 
Lanz & Hubeny (2003, 2007). The table lists the number 
of levels included and, between parentheses, the number 
of super-levels. A super-level groups together several lev- 
els having close excitation energies and assumed to be in 
Boltzmann equilibrium relative to each other. For exam- 
ple, the ion C IV is treated using a total of 25 levels, four 
of them being super-levels. These four super-levels cor- 
respond to the principal quantum number of the external 
electron n = 7,8,9,10, while detailed configurations are 
treated explicitly for n = 2 (s, p) up to n = 6 (s, p, 
d, f, g, h). The abundances are then measured using a 
X^ minimization technique applied to observed spectral 
lines from Tables 2 and 3 and the synthetic spectra. 

In the case of J0720-317 and J2013+400 we also var- 
ied the iron abundance between log (Fe/H) = —6 and —8. 
The effect on the abundance of lighter elements did not 
exceed ±0.02 dex. Therefore, we tabulate only abun- 
dances obtained with models at log (Fc/H) ~ —6. Fi- 
nally, in the case of Feige 24 we fixed the iron abundance 
to log (Fe/H) = -5.5 (Vennes et al. 2000). 

We did not attempt vertical variations of the abun- 
dance of trace element. Chayer et al. (2003, 2006) show 
evidence of a stratification of the oxygen abundance in 
several white dwarfs showing resonance lines of O VI. 
Vennes & Lanz (2001) also found a discrepancy between 
abundance measurements based on O IV and O V in 
Feige 24 and G191 B2B which could be attributed to a 
stratification of oxygen in their atmospheres. 

5.2. Analysis 

Table 8 lists the measured abundances in Feige 24, 
J0720-317, BPM 6502 and J2013+400. New FUSE 
abundance measurements in Feige 24 are consistent with 
the HST STIS measurements of Vennes et al. (2000). 
Figure 6 shows C III line profile fits. The carbon abun- 
dance measurements in J0720— 317 and J2013+400 are 
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TABLE 7 
Model atoms treated explicitly. 







BPM 6502 


Feige 24 


J0720- 


-317, J2013+400 


Element 


I 


II 


III 


IV 


V I II 


III 


IV 


V VI 












H'' 9 


(1) ■•■ 






9(1) ■■• 








He'' 










24 20 








C 




22(5) 


23(7) 


25(4) 


22(5) 


23(7) 


25(4) 




JSfa 




26(8) 


32(7) 


23(8) 




32(7) 


23(8) 


16(6) ■ ■ ■ 


QB. 














39(8) 


40(6) 20(5) 


Si" 




8b 


30(6) 


23(4) 




30(6) 


23(4) 




pb 














14 


17(4) ■ ■ ■ 


S*^ 














15 


12 16 


Fe" 




(36)'^ 


(50) 


(43) (42) 




(43) 


(42) (32) 



'^OSTAR2002, Lanz & Hubeny (2003). 

Earlier version 
''BSTAR2006, Lanz & Hubeny (2007). 



TABLE 8 
Photospheric Abundances for Feige 24, 
AND J2013+400. 



J0720-317, BPM 6502 



Element 




log (X/H) 






Feige 24 


J0720-317 


BPM 6502 


J2013+400 


He 




-3.10 ±0.14 


<-3 


-2.90 ±0.07 


C 


-6.90 ±0.06 


-5.75 ±0.08 


-5.68 ±0.07 


-5.50 ±0.05 


N 


-6.77 ±0.10 


-6.58 ±0.17 


-6.77 ±0.21 


-5.75 ±0.14 


O 


a 


a 




a 


Si 


-6.20 ±0.07 


-6.27 ±0.08 


-6.93 ±0.04 


-6.38 ±0.04 


p 


-7.11 ±0.05 


-7.50 ±0.10 




-8.20 ±0.05 


s 


-6.36 ±0.08 


-6.71 ±0.11 




-6.92 ± 0.05 


Fe 


(-5.5)b 




-7.32 ±0.13 





^See text. 

t'Vennes et al. (2000). 



K 




BPM 6502 



1173 1174 1175 1176 1177 1178 
A(A) 

Fig. 6. — FUSE spectra showing the C III photospheric lines 
of J0720-317, BPM 6502 and J2013+400 compared to the best fit 
models. 

in agreement with the measurements of Vennes et al. 
(1999) based on HST GHRS spectra of the C IV line 
profiles. We could not achieve satisfactory fits of the 



O VI line profiles using the present models. Chayer et al. 
(2003, 2006) explored the problem in a sample of white 
dwarfs showing O VI resonance lines and concluded that 
a reservoir of oxygen concentrated at the top of the at- 
mosphere is able to reproduce the line profiles with plau- 
sible oxygen abundance. We simply note that a similar 
phenomenon is possibly present in Feige 24, J0720— 317, 
and J20 13+400 and we exclude the case of oxygen from 
further consideration. 

Figure 7 shows Si IV and Fe III line profile fits in 
BPM 6502. At Tcff = 20, OOOK and log (Fe/H) = -7.3, 
BPM 6502 seems to follows a trend noted by Vennes et al. 
(2006) suggesting a decreasing abundance with decreas- 
ing temperatures. However, the abundance of iron in 
cooler white dwarfs may depend on their immediate envi- 
ronment. For example, the case of GD 362 at Tcs — 9760 
K and log (Fe/H) = -5.5 (or NLTT44986; Gianninas et 
al. 2004; Kawka & Vennes 2006) may reveal the presence 
of circumstellar debris accreted onto the white dwarf. 
We could similarly argue that the presence of iron (and 
of lighter elements) in BPM 6502 is also a signature of 
accreted material but from the close late-type compan- 
ion. 

Figures 8 summarizes our new abundance mea- 
surements in Feige 24, J0720-317, BPM 6502, 
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Fig. 7. — FUSE spectra showing the Fe III and Si IV (see 
Table 3) photospheric hnes of BPM 6502 compared to the best fit 
modeh 
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Fig. 8. — Measured abundances in Feige 24, J0720-317, 
BPM 6502 and J2013+400 relative to solar values (Asplund et 
al. 2005). The iron abundance in Feige 24 is from Vennes et al. 
(2000). 

and J2013+400. The abundance measurements 

{{X/H)/{X/ H)q) are normalized to solar values (As- 
plund et al. 2005) and presented as a function of the 
atomic number Z. Helium is dominant in J0720— 317 
and J2013+400, but the absence of He I in medium- 
resolution spectra of BPM 6502 obtained by Kawka et 
al. (2000) places a limit on the helium abundance of 
logHe/H < —3. Another interesting contrast emerges 
between Feige 24 and the three other binaries. The case 
of carbon is very instructive (see Figure 6). The abun- 
dance of carbon is « 1% solar in BPM 6502, J0720-317, 
and J2013-I-400, but it is « 0.05% solar in Feige 24. Over- 
all, Feige 24 presents an increased abundance toward 
larger atomic numbers. Although considerable scatter 
is present, one might conclude that, except for P in 
J0720— 317 which appears somewhat more abundant, the 
abundances in BPM 6502, J0720-317, and J2013+400 



are w 0.2-2% solar. Especially in the case of J0720— 317 
and J2013-I-400 this is in stark contrast with the abun- 
dance pattern noted in hot white dwarfs such as Feige 24 
(see Vennes & Lanz 2001) which exhibit a yet unex- 
plained trend showing increasing {X/H)/{X/H)q with 
Z. Although selective radiation pressure, which is very 
effective in the case of elements with a rich line spec- 
trum (e.g., iron), plays an important role in white dwarf 
atmospheres, detailed comparisons between theory and 
observations are not satisfactory (Chayer et al. 1995). In 
the cases of J0720-317 and J2013-f 400, the presence of 
helium (normally not present in hot DA white dwarfs) 
and of a concentration of heavy elements at « 1% solar 
concentration suggest that these elements are accreted 
from the close late-type companion via wind/mass-loss. 
This wind/ accretion scenario was explored in the case 
of the similar white dwarf plus late-type pair 1016—053 
(Vennes et al. 1997). The white dwarf in Feige 24 being 
in a wider binary {P « 4. 2d versus ^ 1 d) and consider- 
ably hotter bears more resemblance to hot, isolated DA 
white dwarfs such as G191-B2B (Vennes & Lanz 2001). 
It is possible to estimate the accretion rate expected 
to bring about the observed level of contamination. 
We adopt the accretion/diffusion model of Michaud & 
Fontaine (1979) as a possible explanation for the abun- 
dance pattern observed in BPM 6502. The model may 
also be tentatively applied to the abundance of light el- 
ements in J0720-317 and J2013-f400. This simple sce- 
nario implies that elements are accreted onto the white 
dwarf surface from the late-type wind; elements heavier 
than hydrogen, the main atmospheric constituent, diffuse 
downward at a certain rate. The equilibrium abundance 
at a certain point in the atmosphere is given by: 

JX/H) ^ M 

{X/H)q M + AirR^pvd 
where (X/H) is the equilibrium abundance at that point, 
{X/H)q is the abundance in the red dwarf mass-loss 
(presumably solar), M is the accretion rate onto the 
white dwarf, R is the white dwarf radius, p is the den- 
sity in the atmosphere at the abundance measurement 
point, and Vd is the diffusion velocity at that same point. 
The white dwarf radius R « 8.4 x lO^cm is given by the 
stellar model, and for the present demonstration we will 
adopt the density p at a Rosseland depth tji = 2/3 — 
p = 7.3 X 10^8 (for BPM 6502) and p = 1.4 x 10"'^ g 
cm-3 (for J0720-317 and J2013-H400) — from the at- 
mospheric models discussed in §4.1. We then calculate 
the diffusion velocity Vd at that location following the 
calculation of Michaud & Fontaine (1979). Table 9 gives 
the diffusion velocity of key elements at r^ = 2/3. 

First, we examine the case of BPM 6502 and estimate 
the accretion rate onto the white dwarf atmosphere us- 
ing the carbon abundance (X/H) = 0.01 x {X/H)q 
and the diffusion velocity from Table 9. A low rate of 
M — 1.1 x lO~^^M0yr~^ is obtained. Next, assuming 
that a constant iron to carbon ratio (solar) is preserved 
in the wind and assuming the same accretion rate, we de- 
termine the predicted equilibrium abundance using the 
correct diffusion velocity for iron. The predicted iron 
abundance is lower than the carbon abundance by the 
factor Vd{C)/vd{Fe) = 0.2. Indeed, the abundance of iron 
in BPM 6502 is about a factor five lower than the car- 
bon abundance. Figure 8 shows that the expected trend 
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TABLE 9 
Diffusion velocities at tjj = 2/3. 



Element 




Vd 

(cm s-i) 




BPM 6502 


J0720-317, J2013+400 


Hell 

cm 

Felll 


0.11 
0.51 


0.003 
0.006 



(X/H)/{X/H)q oc A-^ (full line), where A is the atomic 
weight, seems to apply to the case of BPM 6502. Re- 
newed optical observations of BPM 6502 aimed at placing 
a tighter constraint on the helium abundance would also 
help verify this accretion/diffusion model. The diffusion 
velocity of singly ionized helium is similar to the velocity 
if doubly ionized carbon which would imply similar equi- 
librium abundance. However, the diffusion velocity of 
neutral helium is possibly a factor of 100 larger than for 
singly ionized species which would considerably lower the 
equilibrium abundance of helium (Michaud et al. 1978). 

A similar application of the accretion/diffusion model 
to the cases of J0720-317 and J2013-f-400 implies an 
accretion rate onto both white dwarfs of M — 1.8 x 
lO-i^Moyr-i (He) and M ^ 3.4 x lO^^'^M^yr-^ (C). 
The rates in the cases of J0720-317 and J2013+400 are 
markedly lower than in the case of BPM 6502 possibly 
because of their wider orbital separations or, hypothet- 
ically, their lower late-type mass-loss rates. Vennes et 
al. (1997) determined a similar accretion rate onto the 
white dwarf 1016-053 of M = 8 x IQ''^'^ M^yv''^ based 
on their analysis of the helium abundance in the EUV 
photosphere of that star. 

In summary, the observed abundance pattern in the 
atmosphere of BPM 6502 is consistent with an appli- 
cation of the accretion/diffusion model. The cases of 
J0720— 317 and J2013-I-400 appear somewhat more com- 
plex, but their abundance patterns do bear some similar- 
ities with the pattern observed in BPM 6502. Feige 24 
being in a wider binary than the other binaries studied 
in our sample, the abundance pattern may be dominated 
by ordinary diffusion with possible effects due to radia- 
tive levitation and/or mass loss. In this context, it would 
be instructive to measure the iron abundance in the at- 
mospheres of J0720-317 and J2013-f-400. Such measure- 
ments can be obtained using far UV lines of Fe IV and 
Fe V present at A < 1300 A. 



6. SUMMARY 

We measured the white dwarf radial velocities in the 
close binary systems Feige 24, J0720-317, BPM 6502 
and J2013-I-400 using FUSE spectra. Combined with 
previous optical studies of the binaries, we have up- 
dated the binary properties. Of the four systems, only 
J0720— 317 is likely to come into contact within a Hub- 
ble time and is therefore a representative progenitor of 
the presently observed cataclysmic variables. We have 
also reanalyzed the optical data of Kawka et al. (2002) 
and Morales-Rueda et al. (2005) and showed that the cor- 
rect orbital period of BPM 6502 is 0.3367849 ±0.0000006 
days and that the period reported by Morales-Rueda et 
al. (2005) should be discarded. 

We measured the abundance of trace elements in the 
atmospheres of Feige 24, J0720-317, BPM 6502 and 
J2013+400 and found that the observed abundance pat- 
tern in BPM 6502 can be explained by steady accretion 
at a very low rate (~ 1 x 10^^^ Mq yr^^). In the cases of 
J0720— 317 and J2013+400 accretion from the secondary 
at much lower rates (~ 10^^^ M© yr~^) can explain the 
observed abundances, however the higher abundances of 
P in these stars cannot be explained by accretion from 
the secondary alone. 
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